CNO and pep neutrino spectroscopy in Borexino: Measurement of the deep underground 
production of cosmogenic ^^C in organic liquid scintillator 
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Borexino is an experiment for low energy neutrino spectroscopy at the Gran Sasso underground laboratories. 
It is designed to measure the mono-energetic ^Be solar neutrino flux in real time, via neutrino-electron elastic 
scattering in ultra-pure organic liquid scintillator. Borexino has the potential to also detect neutrinos from the 
pep fusion process and the CNO cycle. For this measurement to be possible, radioactive contamination in the 
detector must be kept extremely low. Once sufficiently clean conditions are met, the main background source is 
**C, produced in reactions induced by the residual cosmic muon flux on *^C. In the process, a free neutron is 
almost always produced. ^^C can be tagged on an event by event basis by looking at the three-fold coincidence 
with the parent muon track and the subsequent neutron capture on protons. This coincidence method has been 
implemented on the Borexino Counting Test Facility data. We report on the first event by event identification 
of in .situ muon induced ^^C in a large underground scintillator detector. We measure a *^C production rate of 
0.130 ± 0.026 (stat) ± 0.014 (syst) day~* ton"*, in agreement with predictions from both experimental studies 
performed with a muon beam on a scintillator target and ab initio estimations based on the **C producing 
nuclear reactions. 

PACS numbers: 25.20.-x: 25.30.Mr: 26.65.+t; 28.20.Gd: 96.50.S-: 96.60.-j 
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I. INTRODUCTION 

Results from solar neutrino (l| and reactor J3l antineutrino 
experiments provide compelling evidence for neutrino oscil- 
lations as the explanation of the long-standing solar neutrino 
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problem |3]. The next goal in solar neutrino physics is prob- 
ing in real time the low energy (< 2 MeV) component of the 
solar neutrino spectrum, which accounts for more than 99% 
of the total flux. This includes neutrinos produced in the pp, 
^Be, and pep nuclear fusion reactions and the CNO-cycle. 
Particularly, pep and CNO neutrinos are an ideal source for 
probing the energy region, between 1 and 3 MeV, at which the 
transition between matter and vacuum dominated oscillations 
is supposed to occur, according to the MSW-LMA oscillation 
solution 1 4]. Furthermore, the pep and pp solar neutrino rates 
are directly related, via the ratio of the cross section of the 
two reactions. Measuring the pep neutrino flux is hence a way 
to study the fundamental pp fusion reaction by which the Sun 
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FIG. 1 : Expected recoil electron energy for different solar neutrinos 
interacting in Borexino assuming 3 year live time exposure, 100 tons 
fiducial volume and a detector energy resolution of 5%/ ^J Eym ev]- 
Neutrino fluxes are derived assuming the Standard Solar Model 
BP2004+LUNA i73|7^ and the LMA oscillation scenario hOil. The 
shaded superimposed area is the expected ^^C background 



burns, and improves our knowledge of the solar neutrino lu- 
minosity, thence yielding a crucial check of the Sun stability 
over a time scale of 10^ — 10^ years by comparison with the 
photon luminosity. CNO neutrinos play a key role on the age 
estimation of the Globular Clusters |5], pivotal in setting a 
lower limit for the age of the universe. 

Deep underground organic liquid scintillator detectors, like 
Borexino and KamLAND, are well positioned to measure pep 
and CNO solar neutrinos. The 1.4 MeV, mono-energetic pep 
neutrinos are particularly well identifiable by the characteris- 
tic Compton-like electron recoil spectrum they produce. The 
main challenge they face is the identification and suppression 
of the ^^C background. ^^C is produced deep underground by 
residual cosmic muons interacting with ^^C atoms in the scin- 
tillator The rate of the process is a function of the location 
and depth of the experiment. As can be seen in fig. ^ the 
^^C background at Gran Sasso falls in the energy region for 
the detection of pep and CNO neutrinos. In 1996, Deutsch |6] 
suggested that ^^C decays could be detected and subtracted 
exploiting the neutron emission in the reaction: 



p (+ secondaries) C — > /i (+ secondaries) C + n. 

(1) 

He proposed using a three-fold coincidence which links the 
parent muon, the neutron capture on protons, and the ^^C de- 
cay. The validity of such technique was studied in detail in 
131 • We apply the three-fold coincidence technique to data 
from the Borexino Counting Test Facility (CTF). This is, to 
the best of our knowledge, the first in situ event by event de- 
tection of ^^C production deep underground. We then use our 
results to evaluate pep solar neutrino detection with Borexino. 



II. " C IN SITU PRODUCTION: THE THREE-FOLD 
COINCIDENCE TECHNIQUE 



^^C /3+-decays with a mean life of 29.4 min and an end- 
point energy of 0.96 MeV: 



"C^" B + e+ + i/e. (2) 



The total energy released in the detector by the decay and the 
following positron annihilation is between 1 .02 and 1 .98 MeV, 
partially covering the best window for the observation of the 
pep+CNO signal (0.8-1.3 MeV). 

The probability to produce ^^C nuclides in muon-induced cas- 
cades was experimentally determined with a target experiment 
(NA54) on a muon beam at CERN 1 10]. The inferred ^^C rate 
for Borexino and CTF is 0.146 ± 0.015 day^^ ton^^ (0.074 
± 0.008 day^^ ton~^ in the pep+CNO neutrino window). 
The study reported in f?'] identified eight different processes 
for the ^^C production in muon showers and provided a quan- 
titative estimate for the rate in all the production channels. 
The result seems robust in view of the fact that the calculated 
production rate matches the rate measured at the NA54 CERN 
facility. 

Two of the production channels identified, ^^C(p,d)^^C and 
^^C(7r+, 7r°H-p)^^C, do not produce a free neutron in the final 
state, and therefore escape any possibility of detection by the 
three-fold coincidence technique. These two production chan- 
nels are referred to as "invisible channels", and they account 
for 5% of the ^^C production rate |7]. 

Neutrons are captured on hydrogen with a capture mean time 
of ^ 250 /is in pseudocumene emitting a characteristic 7 of 
2.2 MeV. Neutrons can also be captured on carbon isotopes 
emitting 7 with larger energy, but the cross section is two or- 
ders of magnitude lower than on hydrogen. 
In order to identify and suppress the ^^C background, each 2.2 
MeV 7 produced in the scintillator from the muon-induced 
showers must be localized in space and time. 
After each muon-induced neutron detection, the three-fold co- 
incidence technique defines a set of potential ^^C candidates 
within a time delay t from the detected muon and inside a 
sphere of radius r from the neutron capture point. We assume 
that no convective currents move the ^^C nuclide from the pro- 
duction point in the time scale of the ^^C mean life. 
In Borexino, the ^^C candidates will be discarded in order 
to increase the pep+CNO signal to background ratio. The 
success of the pep and CNO neutrino measurement will de- 
pend on two main conditions: the minimization of the detec- 
tor mass-time fraction lost to the cuts implementing the three- 
fold coincidence and the achievement of a high efficiency in 
the ^^C suppression. 

The limited size of CTF represents a challenging test for the 
three-fold coincidence technique. The goal in CTF is the mea- 
surement of the ^^C production rate by looking at the time 
profile of the ^^C candidates. 
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FIG. 2: Overview of the CTF detector and of the physical processes 
included in the simulation. 



III. EXPERIMENTAL SETUP 

CTF fs"] is the Borexino prototype detector installed at the 
Gran Sasso underground laboratory. It was designed to test 
the required radiopurity of the Borexino liquid scintillator and 
its purification strategy. The CTF of Borexino was the first 
detector to prove the level of purities needed for solar neu- 
trino physics on a multi ton-scale, in its 1994-95 campaign |9|. 
The active detector consists of 3.73 tons (0.88 ton/m'^ den- 
sity) of the Borexino-like scintillator, a mixture of pseudoc- 
umene (PC, 1,2,4-trimethylbenzene, C6H3(CH3)3) plus 1.5 
g/1 of PPO (2,5-diphenyloxazole, C15H11NO), housed in a 
1 m radius transparent nylon vessel. A 7 m diameter stain- 
less steel open structure supports 100 8" photomultiplier tubes 
(PMT) equipped with light concentrators which provide an 
optical coverage of 21%. 

The detector is housed within a cylindrical tank (11m diame- 
ter and 10 m height) containing 1000 tons of pure water, which 
provides 4.5 m shielding against neutrons from the rock and 
external 7-rays from the rock and from the same PMTs. 16 
upward-looking PMTs mounted on the bottom of the tank veto 
muons by detecting the Cerenkov light in water (muon veto 
system). The veto efficiency is larger than 99.7% for muon 
shower events with energy > 4 MeV. 

A set of analog to digital (ADC) and time to digital (TDC) 
converters records the charge and time information of the 
PMT pulses for each event. During the acquisition, a sec- 
ond identical electronic chain is sensitive to the next event 
occurring within the following 8.3 ms. The electronics can 
therefore detect pairs of fast time-correlated events. The co- 
incidence time between the two chains is measured by means 
of a long range TDC. Further events are ignored until the first 
chain is "re-armed" (^ 20 ms). For longer delays the com- 
puter clock is used providing an accuracy of ^ 50 ms. 
The trigger condition is set by requiring the signal of 6 PMTs 
over threshold within a time window of 30 ns. The corre- 
sponding energy threshold is ^ 20 keV with 50% detection 
efficiency, while 99% efficiency is reached above 90 keV. The 
trigger for the second chain is set at a higher value, corre- 
sponding to 200 keV (99% efficiency). 
The electronic can be also triggered by the so-called after- 
pulses which are spurious pulses following genuine PMT out- 



put pulse. To avoid such effect, the second chain is vetoed for 
20 /is after an event tagged by the muon veto system. 
The energy response of the detector is calibrated run-by-run 
by using the energy spectrum of ^"'C decays, naturally present 
in the scintillator. The measured light yield is ~ 3.6 photo- 
electrons per PMT for 1 MeV electrons. The electronics satu- 
rate at about 6 MeV. 

The position of the interaction vertex is reconstructed by 
means of a maximum likelihood method exploiting the hit 
time distribution. The reconstruction algorithm, calibrated by 
inserting a ^^^Rn source in the active volume, provides a res- 
olution of 10 cm at 1 MeV. 



IV. DATA SELECTION 

The residual cosmic muon flux at Gran Sasso depth (3800 
m.w.e. maximum depth, 3,200 m.w.e. slant depth) has a rate 
of 1.2 m~^ h~^ and an average energy of {E^) ~ 320 ±A-stat 
±llsj,s GeV The requirements in the selection of cos- 
mic muons are two-fold: they must be tagged by the muon 
veto and they must saturate the electronics. Cosmic muons, 
crossing the scintillator, produce enough light to blind the de- 
tector 

For each detected muon, we select the following event in the 
time window T„ = [20, 2000] /is as a candidate event for a 
neutron capture 7. The probability that a random event (R = 
0.04 s^^ rate) is detected instead of the 2.2 MeV 7 has an up- 
per limit equal to T„ x R ~ 8 x 10^^. We measured the mean 
capture time of neutrons on protons equal to 257 ± 27 /is, tak- 
ing into account also events with double neutron emission. 
For each muon-gamma coincidence, ^^C candidates are se- 
lected in a subsequent time window T^ = 300 min, 10 times 
the ^^C mean life. 

Random coincidences collected in this window are mainly 
2iOBi (Q;3=l-16 MeV) and ^^K (Q/5=1.32 MeV BR=0.893 
and Q£;c=l-51 MeV BR=0.107) contamination and external 
7 radiation, while ^^''Bi (Q/3=3.27 MeV) events are discarded 
through the ^^^Bi-Po coincidence. 

The time profile of the background is expected to be flat on the 
scale of 300 minute since the background rate is constant and 
random coincidences are not correlated with cosmic muons. 
The only bias is introduced by the end of the data run (typ- 
ically lasting 2-3 days) which interrupts 8% of the selection 
windows. In such cases the window is completed to 300 min 
from a random instant in the run in order to correctly maintain 
the background time profile flat. We estimated that the corre- 
spondent probability to loose a ^^C event is 1% fT^ . 
The definition of the optimal energy range of observation, 
1.15-2.25 MeV, to detect the ^^C decays, depends on two main 
requirements: the enhancement of the signal (^^C decays) to 
background (random events) ratio and the minimization of the 
systematic errors introduced by the energy scale uncertainty. 
In case 7's from the positron annihilation escape the vessel 
and deposit energy in the water buffer, the detected energy of 
the ^^C decay falls below the observation range. Defining a 
0.8 m radius fiducial volume, we reduce non-contained events 
by a factor 20. Further, the radial cut avoids distorting optical 
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FIG. 3: Scatter plot of the simulated C-decay radial position 
vs energy. Solid lines represent the cuts applied in the analysis. 
Events with an energy lower than 1.022 MeV are due to the escaping 
positron annihilation 7's. 



effects on the border like the total reflection due to the differ- 
ent refractive indexes of the scintillator and the buffer 
The last applied cut exploits the spatial coiTelation between 
the ^^C and the neutron capture points. The events are in fact 
selected in a sphere of radius r centered on the reconstructed 
2.2MeV'-f: for r = 35 cm the background is suppressed by a 
factor larger than 20 while the signal is reduced only by a fac- 
tor ~ 2. 

The efficiencies and optimal parameters of the cuts here dis- 
cussed have been quoted via the Monte Carlo simulation de- 
scribed in the next section. 



the neutron capture on hydrogen. 

The CTF code simulates in detail the detector geometry 
including the nylon vessel and the phototubes. Each en- 
ergy deposit is converted into optical photons which are 
propagated inside the detector until they are absorbed in the 
detector material or detected on the PMT's. 
The tracking code provides a detailed simulation of the main 
optical processes like the scintillation light production, the 
absorption and reemission processes in the scintillator and 
diffusion on the nylon vessel. 

After all, the same reconstruction code used in the real data 
introduces the energy and spatial resolution effects on the 
simulated ones. The final ^^C radial and energy spectra are 
shown in figure|3 

The simulated neutron capture mean time, 254 ± 1 /is, is in 
good agreement with the measured one, 257 ± 27 /is. 
Secondary particles generating a ^^C event without trig- 
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V. THE MONTE CARLO SIMULATION 



FIG. 4: Cumulative probability to contain a 
radius r centered on the 2.2MeV'Y 



C event in a sphere of 



An accurate quantification of the cut efficiencies requires 
a full simulation of the ^^C production process from the 
muon-induced showers originated in the rock to the neutron 
capture and to the ^^C decay. 

The Monte Carlo has been developed in two main steps. 
First we generated and tracked muons and the subsequent 
cascades with a FLUKA-based code 1 13]. The code simulates 
a 320 GeV muon-beam, downward oriented and uniformly 
distributed over the entire CTF water tank. At this step, the 
geometry is simplified to only four volumes: 4 m of rocks 
(CaCOa and MgCOa) 1 14|, the air, the water of the CTF tank 
and finally the scintillator as shown in Figure|2] The purpose 
of the FLUKA-based simulation code is the generation of 
neutrons in scintillator and their propagation in the whole 
detector. 

In the second step, an ad hoc code, named CTF code 
generates, tracks and reconstructs ^^C decays and 
2.2MeT^7's from the neutron capture. The coordinates of 
the neutron production {Pp) and capture (Pc) points from 
the FLUKA output are input parameters in the CTF code: 
Pp corresponds to the origin of the ^^C-decay while Pc is 
assumed as the starting position of the 2.2Af eVj produced in 



gering the muon veto have been investigated. From the 
simulation, we expect mainly 7's (91.8%) and e+-e~ pairs 
(8.1%). Their contribution to the invisible ^^C production 
rate has been estimated in less than 5x10^^ day^^ (99.99% 
C.L.) by convoluting their rates with the ^^C production cross 
sections |7]. 

The main inefficiency in the measurement is due to neutrons 
escaping the vessel. If the neutron, indeed, is captured in 
water and the subsequent 7 does not deposit energy in scintil- 
lator, the /i — 72.2A/ey coincidence is not triggered and the 
signal is lost. Neutrons escaping the 1 m CTF vessel account 
for 26.8%. For ^ 50% of the fully contained neutrons, the 
associated ^^C event falls in a 35 cm radius sphere centered 
on the reconstructed 2.2A/eI^7, as shown in Figure|4] 
All the cut efficiencies are quoted in TableHl 



VL THE DATA ANALYSIS 

The analyzed data set corresponds to an effective detector 
live time of 611 days (June 2002, February 2005). 
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The time profile of the data sample selected by the three-fold 




Coincidence Time [min] 

FIG. 5: Fit of the data sample time profile selected by the three-fold 
coincidence (x^/ndf = 9.7/13) with free parameters A and b (see Eq 

m 

coincidence technique, shown in Figure|5] is fitted with: 

Pit) = -e-^ +b, (3) 
r 

where the free variables in the fit, A and t, are the number 
of ^^C nuclides and the ^^C mean life, respectively. The fit 
finds r = 27 ± 11 min (A = 53 ± 13, bxT^ = 166 ± 17 
and x^/d.o.f = 9.7/12) in agreement with the nominal value 
(29.4 min), proving the robustness of the three-fold coinci- 
dence technique. Moreover, if the 300 min window is started 
independently from the /i — 72.2A/ey coincidence, the fit is un- 
able to identify any feature compatible with a decay function. 
Performing the fit with t fixed to the nominal value, the ^^C 

TABLE I: Efficiencies for the ^^C production rate measurement in 
CTF. 



Efficiency Reason 


Value 




visible channels 


0.955 


£ end 


end of run 


0.990 


£t 


/i-2.2A/eV7 coincidence 


0.925 


£ escape 


contained neutrons 


0.732 


Ec 


'^C energy E e [1.15, 2.25] MeV 
2.2MeV7 energy E > 0.2 MeV 
"C-2.2Me\/7 distance d < 0.35 m 


0.563 


Total 




0.360 



production rate is computed from: 

= ^4—-; r 

2 TTT pT £vis ' ^end ' ' ^ escape ' 

(4) 

= 0.130 ± 0.026{stat) ± Q.014:{syst)day-^ton-^ 

(A = 54± ll,bxT„= 164± 15andx^/d.o.f=9.7/13)where 
r is the selected volume radius (0.8 m), p the scintillator 



density (0.88 g/cm^) and T the detector live time (611 days). 
All the efficiencies in Eql^are reported in TableU 
The systematic error has been derived by propagating the 
uncertainties of the reconstruction position (^ 1.5%) and of 
the light yield (^ 8.5%) in Eq0] The systematics takes also 
into account the stability of the result when the cut parameters 
vary around the optimal values. 

The analysis measured rate is in good agreement with the 
expected one from the CERN experiment: 0.146 ± 0.015 
day^^ ton^^. 



VII. DISCUSSION 

The success of the three-fold coincidence technique in se- 
lecting ^^C events and in evaluating correctly their production 
rate is promising in prospective of deep underground liquid 
scintillator detectors. 

The expected rates for pep and CNO neutrinos in Borex- 
ino are 0.021 and 0.035 day-^ ton^^ (BP2004h-LMAh-LUNA 
E US 113]), respectively. In the energy range of observa- 
tion [0.8,1.3] IVIeV, beyond the ^Be-i/ electron recoil energy 
spectrum, the pep+CNO signal. Si,, is reduced to 0.015 day^^ 
ton^^. In the same window, the expected contamination from 
"C is then about 5 times higher (Bn^ = 0.074 ± 0.008 day-^ 
ton^i). 

A second background contribution arises from the trace con- 
taminants in the scintillator mixture. Assuming for the ^ssy 
and ^'^^Th a concentration level of lO^^^g/g and lO^^^g/g for 
the "'^*K, the non-cosmogenic contaminants, B„ c., contribute 
to the pep+CNO window with 0.006 day^^ ton^^ 
In order to reach a signal-to-background ratio equal to 1, 
the detection efficiency of Borexino must be larger than 1 - 
S,y/(Biic'H-B„.c.) = 0.81. The detection efficiency is limited by 
the physics (^^C(X,Y)^^C invisible channels), by the detector 
itself (low energy threshold and dead time between sequential 
triggers) and by the software cuts in time and space around 
the neutron capture 7's. Since, in fact, the three fold coinci- 
dence does not identify the single ^^C decay but localizes it 
in a spherical volume Vii(^, the entire volume Vnc must be 
discarded for a time equivalent to few ^^C lifetimes. Thus, the 
main challenge will be the minimization of the detector mass- 
time fraction loss. 

Assuming a neutron rate of 1.5x10^^ m^^ f?! I2TI1 . we 
estimate that, even including the trace contamination, Borex- 
ino can reach a signal-to-background ratio equals to 1, loosing 
only 14% of the data f7l l2^ . The optimal cuts and the relative 
efficiencies expected for Borexino are quoted in Table IIII 
Furthermore, the Borexino collaboration is investigating i23ll 
the possibility to improve the three-fold coincidence tech- 
nique by exploiting the muon track: the reconstruction of the 
muon track leads, in fact, to the definition of a cylindrical vol- 
ume around the track itself. Intersecting the cylindrical vol- 
ume with the spherical one centered on the 2.2MeV'y, Borex- 
ino can efficiently remove ^^C events while reducing signifi- 
cantly the fraction of data loss. 
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TABLE II: Predicted ^^C decay detection efficiencies for the Borex- 
ino detector in order to reach a signal (pep+CNO vi) to background 
(^^C decays) ratio equal to 1. The trace contaminanion is assumed 
at the level of 10""g/g for ^^'^^U and ^^^Th and 10~^'^g/g for ""'K. 



Efficiency Reason 


Value 


e-uis visible channels 


0.955 


£( jj.-2.2AI eV"f coincidence 


0.989 


£c 2.2MeV'y energy E > 0.2 MeV 


0.954 


£a ''C-2.2MeV'y distance d < 1 m 


0.984 


Et ''''C-2.2MeVj coincidence time T < Sxrii^ 


0.993 


Total 


0.880 



dence technique is a powerful tool for isolating and discrimi- 
nating the ^^C background. 

The results also indicate an agreement with the theoretical cal- 
culation in |7]. When combined with the prediction that the 
overall rate of ^^C produced without free neutrons in the final 
state is limited at 4.5%, this observation indicated that Borex- 
ino should be able to minimize the ^^C background at a level 
compatible with the observation of pep neutrinos. 
In prospective of Borexino, such result opens a new window 
in pep and CNO neutrino spectroscopy. 
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get experiment [10], demonstrated that the three-fold coinci- 
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